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Abstract 

We give a complete treatment of factorisation of Deeply Virtual Compton Scat- 
tering (DVCS) in the generalised Bjorken limit, using the local Operator Product 
Expansion (OPE). The method allows a straightforward proof that, at leading twist, 
the DVCS amplitude factorises into an integral over coefficient functions and Skewed 
Parton Distribution Functions (SPDFs). The integral is well defined for on-shell fi- 
nal state photon if the Wilson coefficients satisfy a certain factorisation condition, 
which we derive. We also show that it enables a simple proof that soft singularities 
either cancel out or, in the case where the final state photon is on shell, are inte- 
grable. This confirms the argument of Collins and Freund. Further, we repeat the 
tree- level calculation of twist-three contributions to DVCS off a scalar target, where 
factorisation was found to be violated. We propose a new definition of the structure 
functions and calculate the coefficient functions, which are such that factorisation 
works. 
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1 Introduction 



In recent years, there has been much interest in the process known as Deeply Virtual 
Compton Scattering (DVCS) [0, 7*p — > jp. In the generalised Bjorken regime, where 
the invariant mass of the virtual photon 7* is very large compared to all other kinematic 
scales except the 7p centre-of-mass energy squared, it has been shown 0, |, ^, |[ that the 
amplitude for DVCS factorises into a perturbatively calculable coefficient function and 
a non-perturbative factor, the Skewed Parton Distribution Function (SPDF) 0. This 
factorisation generalises that of the forward Compton amplitude relevant to Deep Inelastic 
Scattering (DIS). Thus, measurements of DVCS can be used to extract information on 
SPDFs. One-loop perturbative calculations have been performed P, |], ^ for some of the 
possible Lorentz structures which demonstrate factorisation. 

Our principal purpose will be to provide a complete treatment of the factorisation of 
the leading term for the DVCS amplitude from the point of view of the local Operator 
Product Expansion (OPE). Previous proofs of factorisation have applied rather different 
methods. In refs. 0, factorisation was considered by analysing the mass singularities 
that appear when all small scales in the problem go to zero. All collinear mass singularities 
can then be absorbed by the SDPFs. However, soft singularities, where soft lines connect 
the hard coefficient to the SPDF, were shown to make leading order contributions. For 
factorisation not to be spoiled one would have to show that the soft singularities are in 
fact integrable This was shown explicitly to be the case at the one-loop level, given 
that the SPDFs are sufficiently smooth. For the all-order proof |^ , a general argument 
can be applied. An important feature of the local OPE treatment of the present paper is 
that one can show quite easily that all such soft singularities either cancel out or, where we 
put the final-state photon on shell, are integrable, given the same smoothness assumption 
for the SPDF. Moreover, the local OPE provides a concise and elegant scheme in which 
to perform calculations. 

More recently, the effect of contributions from twist-3 SDPFs has been considered 



11] , |1^]. It has been shown that, even at tree-level, they are such that the factorisation 
formula diverges. The second purpose of this article is to revisit these calculations. We 
shall find that a suitable redefinition of the structure function to which the twist-3 SPDFs 
contribute will be sufficient to give them a factorisation formula which does not diverge. 

In order to demonstrate how factorisation works in the local OPE method, we will first 
consider an analogue of DVCS in scalar field theory. This will enable a demonstration of 
the method without a plethora of indices. The local OPE analysis is a straightforward 
extension of that applied to DIS. In DVCS we must include total derivative operators 
in the OPE, which have non-vanishing off-forward matrix elements since the final state 
proton possesses non-zero momentum relative to its initial state. Furthermore, the local 
OPE is an expansion in an unphysical region, as in DIS. We will use a dispersion relation 
to reconstruct the full, factorised leading twist amphtude. We shall find that dispersion 
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relation integral is only defined provided that the Wilson coefficients obey a factorisation 
condition. 



We then move on to consider the local OPE for the Compton amplitude at tree level 
in QCD. The organisation of the OPE into terms each obeying current conservation is 
demonstrated. In particular, for the case of an unpolarised proton, we propose a suitable 
basis of four structure functions, Tt, Tl and T^, each multiplying a Lorentz structure 
which respects current conservation. Since the OPE is an expansion of two conserved 
currents, we shall check that the OPE is indeed consistent with current conservation. 
However, as recently discussed in ref. [|T0|, the current conservation is only satisfied be- 
cause of a cancellation of terms between parity-even and -odd operators. Fortunately, 
it turns out that the twist-three operators can, at tree level, be related to the twist-2 
ones |]lOl- Thus, the OPE can be reorganised into two pieces of even and odd parity, 
each respecting current conservation. The twist-3 pieces contribute exclusively to the 

structure functions. We calculate explicitly their coefficient functions and show that 
factorisation is not violated. This is because the corresponding Wilson coefficients obey 
the factorisation condition mentioned above. 

The article's contents are as follows: Section ^ briefly describes the choices of kinematic 
variables used for DVCS; section |^ contains the treatment of the local OPE and dispersion 
relation used to reconstruct the factorised amplitude from the local OPE, together with 
a simple model for the Wilson coefficients from which are derived the corresponding 
coefficient functions; in section ^ we consider the OPE at tree level in QCD and find the 
tree-level predictions for the four structure functions in terms of SPDFs; finally we give a 
summary in section ^. An appendix on SPDFs is provided at the end. 



2 Kinematics 



In this section we give some definitions of kinematic variables. In particular, we define 
the generalised Bjorken limit in which the amplitude can be shown to factorise. Referring 
to fig. ^ we label the initial and final state photons' momenta by qi and q2 respectively 
and likewise the initial and final state proton's momenta by pi and p2- In terms of these, 
we define three independent 4-vectors: 

P = \ (Pi +P2), (1) 

A = P2-Pl- 

The convenient, Lorentz- invariant kinematic variables are as follows: 
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Figure 1: Definitions of momenta for DVCS. 



t 

X 



2p ■ q' 
q ■ A 
2q ■ p 



(2) 



The generalised Bjorken limit is defined by taking large, while keeping x and ^ 
finite. In this limit, we can neglect any quantities suppressed by powers of t/Q"^ and 
M'^/Q'^, where pl=pl = M'^. 

For future reference, we also write down the s- and u-channel variables together with 
the invariant masses squared of the photons: 

s = {q + p) = + - -t, 

^ X 4 

u = {q-pY = --^ ^ + - -t, 

^ ^ X 4 

2 QHx-n 1 

ql = - ^^^^ + 4 t, (3) 

Q' + ^ 1 

92 = + T ^• 

X 4 

Since Q'l = for DVCS, then we have a constraint on the four, otherwise independent 
variables in eqs. (0): 

However, for the purposes of the local OPE analysis in section |^, we relax this constraint. 
The photon always can be put on shell (g| = 0) at the end. Moreover, for other choices 
of g| one can describe processes other than DVCS. For example, if g| > 0, we can also 
describe, for example, dilepton production 7*p —>■ '~f*p —>■ l^l~p. 

In the next section we show that, in the generalised Bjorken limit, the amplitude 
factorises, to leading order in l/Q^, into a convolution of a perturbative coefficient with 
a SPDF. 
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Figure 2: Branch cuts in the complex x' plane. 

3 Construction of Compton Amplitude from the Lo- 
cal OPE 

We now give the proof of factorisation of the leading twist contributions to the Compton 
amplitude in terms of the local OPE. A brief outline is as follows: As with the local 
OPE treatment of the total DIS cross-section, where ^ = t = 0, the local OPE gives 
a simultaneous expansion in inverse powers of and inverse powers of x. The series 
thus converges only in the unphysical region |a;| > 1 where the Compton amplitude has 
no discontinuities (physical cuts). However, one can use the analytic property of the 
amplitude to reconstruct the full amplitude amplitude via a dispersion relation in x. The 
amplitude T can thus be expressed, to leading order in by the following convolution 
formula 

T(x, e, Q\ t) = C - c{-, ^, Q') f{u, e, t), (5) 

J-1 u u u 

where C are perturbatively calculable coefficients which can be derived from the leading 
Wilson coefficients in the local OPE and where / are the SPDFs. This factorisation inte- 
gral is finite for DVCS, ^ = —x, provided that the Wilson coefficients satisfy a particular 
condition which we state explicitly below. This condition is of particular importance since 
factorisation does not work if the Wilson coefficients do not satisfy it. When we come to 
consider twist-3 contributions in section ||, we shall need to define the structure functions 
in such a way that it is satisfied. 

In one of the papers which proved the factorisation formula (|]), Radyushkin showed 
that the coefficient functions are free from soft singularities. In this context, let us take 
soft singularities to be any divergences in the coefficient function C as u ^ ie. at 
the endpoints of the ERBL region of the SDPF (see Appendix) where either parton line 
in fig. ^ becomes soft. Such soft singularities would take the form of terms proportional 
to 1/(m^ — ^^). In the limit u — > ±.^, all of these should cancel out. We shall show this 
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Figure 3: Factorisation of DVCS amplitude in the generalised Bjorken limit. The partons' 
longitudinal momentum fractions are labelled. 

below explicitly, starting from the local OPE. The literature has also been concerned 
with additional soft singularities which should appear as the final-state photon goes on 
shell, ^ = —X. We will show that these are due to terms proportional to l/{u'^ — x'^) 
which tend towards — C^) as a; — *• It was also shown [Q that the singularities 

are integrable and therefore pose no threat to factorisation, given that the SPDFs are 
sufficiently smooth. (See Appendix.) The local OPE enables us to verify that this is case. 

At the end of this section, we consider a simple model for the coefficient function to 
show explicitly the cancellation of soft singularities in the Compton amplitude. This will 
enable the reader to follow the local OPE formalism in a simple example. 

To start, we consider the following amplitude in scalar field theory (fig. |l]): 



T(g,P,A) 



iq-z 



P2 



T(02(iz)02(-K 



Pi 



(6) 



The amplitude T(Q^ t, x, 

is an analytic function of its arguments which contains branch 
cuts in the physical regions. Since we are going to require later a dispersion relation in 
the variable x, we are interested in the branch cuts in the complex x plane which appear 
where the following invariant masses are positive: 



s > 
u>0 
ql>0 
ql>0 



< X < 1, 
-1 < a; < 0, 

(0 < X < OR (^ < X < 0), 
(0 < X < -0 OR (-^ < X < 0). 



(7) 



The inequalities are neglecting the corrections of 0(M^/Q^) and 0{t/Q'^) from eqs. 
These regions are shown as two wavy lines in the x-plane in fig. 0. 



5 



The dispersion relation is derived by using Cauchy's theorem in a region away from 
any of the abovementioned discontinuities in the x plane (fig. |^), we have 

e, Q^ t) = -P^ T{x\ e, g^ t). (s) 

Deforming the contour 7 to 7', and under the assumption that the contour at infinity 
vanish we have 

T(x,e,Q^t) = -- / ^^PTix',^,Q^t), (9) 
vr J-i X — X 

where 

2^pT(x^e,Q^^) = -discT(x',e,Q',t). (10) 

x' 

The 'discontinuity' is defined as disc/(z) = lime^o(/(^ + ^e) ~ f{z — ie)). One can see 
from fig. that the discontinuity includes contributions from s- and u-channels and also 
from the channels associated with the invariant masses of the photons, qf^2- The latter of 
these of course vanish in the forward (DIS) limit ^ ^ 0. 

Since we shall below be obtaining an expansion for T at large x from the local OPE, 
we expand the dispersion relation for large x: 

Thus, in an expansion at large we can identify the coefficients with jth. Mellin 
moments of the spectral function pT- So, the steps required in constructing the factorised 
amplitude are 

1. Use the local OPE to obtain the large x expansion coefficients and identify these 
with the Mellin moments of the spectral function. 

2. Invert the Mellin moments to obtain the spectral function. 

3. Reconstruct the full, analytic amplitude using the dispersion relation (^. 

To implement step one, let us now use the local OPE to obtain the expansion at large 
I a; I and large Q^, keeping the leading terms in 1/Q^. We do this for the simple case of 
scalar fields, to show the general method without any indices. The time-ordered product 
of the currents in (|^) has the following expansion: 

T(02(i;.)02(-i^)) ^ (12) 
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The Wilson coefficients C™'-'(z^) are perturbatively calculable and will give rise to the 
'hard' or Q^-dependent part of the factorised amplitude. The neglected terms are higher 
order in and are small as 0. In momentum space, this limit corresponds to 

oo. Notice also the presence oi m < j total derivatives acting on the spin j — m 
local composite operators. If we take off-forward matrix elements of these operators, we 
can write them in terms of reduced matrix elements as follows: 



P2 



i d ■■■ i d^'^: (0) 



Pi) = J2 lA^i ■■■ iA^»p'^'+^---p'^^ (O^'-'). (13) 



even i=0 



By crossing symmetry, only even powers of A are included in this decomposition. Taking 
the Fourier transform of the off-forward matrix element of the expansion (|l^) and using 
(|T3]) we get the amplitude in momentum space 

T(x,e,Q^t) "'^"^ y: x-^ c^HQ')r E e{o'''-n, (14) 

even j=2 even m=0 even i=0 

where 



d_ 

Comparing eqs. (14) and (|lT]) we can read off coefficients of x^^ and we get 

"1 1 ^ _ 

even m=0 even i=0 



C™'i(g2) = (q'^) j d^ze'"-' C^^^{z^). (15) 



- dxx^-' PT{x,i,Q\t) = - E C^'KQ^)r E C{0^'''n- (16) 
vr Jo / — n — n 



Eq. (0) expresses the factorisation of the moments of the spectral function. The hard, 
perturbative physics is encoded in the Wilson coefficients, (7, while the non-perturbative 
physics is in the reduced matrix elements. Notice that the reduced matrix elements appear 
in polynomials of ^ of order j — m. In the appendix on SPDFs it is shown that these 
polynomials are equal to moments of the SPDFs: 



fiU) = f'duu^-' f{u,^,t) = I t: e{0'''). (17) 

Jo I „,.Z~-r. 



- f'dxx^-' PT{x,^,Q\t) = E C^HQ')C f-'^iU)- (18) 
71 Jo rrl-n 



even i=0 

With this identification, (^) becomes 

dx x^ ^ PT{x,^,Q'^,t) = 

even m=0 



For step two, we need to invert the Mellin transform in eq. (0). One can verify that 
this is solved by the following expression: 

-pT{x,^,Q',t) = 9{x>0) f'-C{-J,Q') f{u,U) 
vr Jx u u u 
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- e{x<o) r -Ci-A,Q')fiu,^,t) 

J-1 u u u 

- \d{Q<x<\i\) f - D{-, i,Q2) f{u,^,t) 
2 J-i u u u 

+ 7J^(-Iei<^<0) / -/^(-,-^,Q^)/(n,e,t), (19) 
2 j-i u u u 



where the coefficient functions C and D have been introduced so that 

»i 







x^-^ Cix,^,Q' 



even m=0 
oo 



„1 oo _ 

Jo _,._r~:_a 



(20) 
(21) 



even m=2 



Notice that we have four terms in the spectral function. These correspond to the s-channel, 
u-channel and g^/2"Channel cuts respectively. 

Finally, for step three, we can insert (0) into the dispersion relation (P) to obtain our 
final form for the factorised amplitude: 



T{Q',t,x,0 



where 



c'(Q^ _ , 



X 

u' u' 



2 Jo 



1 du 

-1 u 



dx' 



5 5 

u u 



(22) 



1 1 
+ 



- / dx' 
2 Jo 



1 



x' + X 

1 



u u 



x' — X 



x' + X 



u u 



(23) 



Three points are now in order: First, the soft singularities can be shown to cancel out 
for X 7^ ±^: The key point about the infinite series in (pO|)-(|2T|) is that they diverge as 
|.^| — > 1 so that they each behave as 1/(1 — S,'^). In the convolution formula in eq. ([19D , 
these singularities appear as m ^ so that, typically. 



u u 



u u 



(24) 



Since the ERBL region of the SPDFs is— |^|<m<|.^|,C and D diverge at the endpoint 
of the ERBL region. However, in this case the singularities cancel between the functions 
C and D. We can check this by taking the limit u — of the coefficient function C to 
get 



hm C{-J,Q') 



— / dx' 



X' — X 



+ 



x' + X 



limfc(^,^,Q2) 



.(25) 
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C and D both have the same divergent behaviour, as can be deduced from the infinite 
series fEOD-fCT), which cancels between them. 



The second point concerns the presence of terms that go hke — x^) which arise 

from the first integral in eq. (0). For example, if C{x,C,) = 5(1 ~ x), we would have 

^(-,-) = -^^- (26) 
u u — x-^ 

These are the additional soft singularities which, as mentioned before, arise as either 
photon goes on shell, ^ —>■ ±x. It is quite easy to show that we will not get terms any 
more singular than in eq. (|26|) . As mentioned above, in the DVCS limit x —>■ — ^, this kind 
of soft singularity approaches the endpoints of the ERBL region u = However, given 
a smooth enough SPDF (see appendix), the singularity remains integrable. We could not, 
however, tolerate terms more singular than in eq. (p6|). 



The third point is that if we put one of the photons on shell by taking the limit 
^ — *• ±x in the factorised amplitude (P^D, it appears that the coefficient function might 
diverge logarithmically. This happens when either limit of the second integral over D 
in eq. ( PBD is trapped at either of the two poles in x ± x' and indicates the presence of 
ln^(l ±^/x) terms. However, no such difficulty occurs provided that D{x,C,,Q^) ^ as 
^ — *• ±x. A sufficient condition for this to happen can be derived in terms of the Mellin 
moment of D in eq. (pT]): As j becomes larger, the integral samples the region = ±x 
with a higher weighting. If D{x, ^, Q"^) goes like ^ — x in this region, one can show that 

dx x^-^ D{x^,^,Q^) ~ dx x^'^ {I - x) 
Jo Jo 

(27) 



1 J— *oo 1 



j(j + 1) 3 

Thus, from the r.h.s. of eq. (pT]), a sufficient condition for the absence of ln^(l ± ^/x) 
terms is 

C^+3,j 3^ j__ (28) 

Eq. (^) is crucially important to factorisation as treated in the local OPE: If it is 
satisfied then factorisation is valid in the DVCS kinematics. All NLO calculations [0, |^, Q 
of twist-2 contributions that have been presented are of course consistent with (|28|). 



A simple model 

To illustrate the application of the local OPE, let us now consider a simple example which 
has all the features of a factorisable amplitude. We choose a simple but non-trivial model 
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for the Wilson coefficients which satisfy eq. (^81) . Then we perform steps two and three 
in the method described above to obtain the factorised form for the 'amphtude'. 

Suppose, then, that we performed an OPE and obtained the following expansion 

T(a;,0 ^^"^ x-^ E r C""'^' f -'"(O, (29) 

even j=2 even m=0 

where 

C"^'^ = -, (30) 

and f-'iC) is the jth. moment of the skewed PDF. 

For step two, we take the coefficients of and perform the inverse Mellin transform 
w.r.t. j to obtain: 

1 fill T f 

-Pt{x,0 = 0{x>0) / -C{-,^)f{u,0 

TT Jx U U U 

pi rjii qr f 

- e{i>x>Q) / _/}(-,!) + ■■■ , (31) 
Jo u u u 

where the neglected terms are those which are required to make pt an odd function of x. 
We also take ^ > 0. The coefficient functions are given by0 



oo 



C{x,i) = E ra:(l-x) 

even m=0 

= u{l-u)j^^, (32) 

even m=2 ^ \ ^ / 



1 - 7 1 7^- (33) 



For step three, we replace pr into the dispersion relation. This gives us a double 
integral. Performing the integral over x' we arrive at the final convolution integral: 



TiM) = -\f- 

2 J-i u 



' X , (l — u/x\ x'^ , u^W 
1 + - In . . ' - — In 1 - — * * 



u \l + u/x) \ x"^ ) ) — ^"^ 



i \l+^/x/ \ x"^ ) ) u 



2 



(34) 



The coefficient of / in the integral can be seen to have the property that it is free of soft 
singularities, ie. as m — >■ It is also finite in the DVCS kinematics ^ = ±x. 

^The inverse Mellin transform of l/((j + !){] + 2)) is the function u{l — u). 
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4 The Local OPE in QCD 



In the previous section, we applied the local OPE to the analogue of DVCS in scalar 
field theory. We now consider the QCD case. First, it will be important to make some 
comments about how to organise the local OPE. It is important to keep conservation 
of both electromagnetic currents explicit. To this end, we shall start the section by 
writing down a general decomposition of the Compton amplitude in terms of Lorentz 
scalar structure functions, keeping the photons of unequal mass. For simplicity we shall 
retrict our attention to the case of a scalar target which, for practical purposes, can 
be an unpolarised proton. We also work in the approximation of neglecting any terms 
suppressed by 0{t/Q^). All structure functions are such that electromagnetic current 
conservation is obeyed for each. 

Afterwards, we shall derive the local OPE for the product of electromagnetic currents 
at the tree-level approximation. We show explicitly that current conservation is satisfied 
because of a cancellation between parity-even and -odd, twist-3 operators. This has been 



emphasised in ref. ||T0[- The first purpose of this section is to review briefiy their work. 
Another important observation of ref. was that the operator identities which allowed 
the twist-3 terms to cancel also enable us to relate the twist-3 operators to purely twist-2 
terms. Once the OPE is recast in terms of these twist-2 pieces, it can be seen that even 
and odd parity sectors satisfy current conservation separately. This means that the OPE 
does indeed enable us to partition the operators of opposite parity from each other in 
different Lorentz structures. 

Once the OPE has been treated, we will then take off-forward matrix elements of 
the OPE expansion in unpolarised proton states. This will enable us to give tree-level 
predictions for the structure functions in terms of integral convolutions of twist-2 SPDFs 



with coefficient functions. This has also been achieved in the literature before pQ| , [11], 
However, the structure functions which receive contributions only from twist-3 pieces were 
found to diverge logarithmically in the DVCS limit. In our treatment, we give a different 
definition of the structure functions for which the tree-level convolution integrals are finite. 



Structure functions 

Let us write down the decomposition of the Compton amplitude tensor T'^'^ for the case 
of a scalar target. We only include those Lorentz structures which are of even parity and 
satisfy current conservation, ie. 

{q+kA),T''%q,p,A) = (g-iA).T'^'^(g,p,A) = 0. (35) 

It will be convenient to introduce the following four-vectors: 

q^' = qi" + xp^, (36) 
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(37) 



These are such that q ■ A_|_ = while (f & A^^ = 0{t). Thus, keeping only terms of order 
Ax and neglecting orders A^ and t, the most general structure function decomposition 
is as follows :0 



T^'^(g,p,A) 



— ip^^q"^ + 2xp^p'' + |/'A^)Wr(x,e,Q',t) 
I ■ p \ ' j 



TL{x,^,Q\t) 




2q ■ p 
0{Al,t) 



T3-(x,e,g^^) 

T3+(x,e,Q^^) 

(3^ 



The structure functions and are both even under crossing symmetry and are 
therefore even functions of x and ^. Tt and T-^ are odd under crossing symmetry and are 
odd functions of x and even functions of ^. In the forward (DIS) limit, ^ & Aj^ — >■ 0, the 
above decomposition reduces to just two terms with Tt and Tl, the familiar transverse 
and longitudinal structure functions. 

The reader may object however that, when we look at the case of DVCS, ^ = —x or 
equivalently ^ = x, the Lorentz structures multiplying Tl and T^ become large. Other 
authors have not used these definitions for this reason. However, the divergence is physical 
and should not be included in the definition of the structure functions. As a result, the 
factorisation formulae previously found 0] for the alternative definitions of T^ have been 
found to be divergent. The divergence does not appear in physical amplitudes because, 
at this order in A^, the physical polarisation vectors of the final-state real photon are 
perpendicular to the vector 



p" + 



x + i 



(39) 



Local OPE at Tree Level 

In order to obtain predictions for the four structure functions in terms of skewed PDFs to 
leading order in the strong coupling, we now use the OPE method introduced in the last 
section to demonstrate factorisation. What we shall find is that vanishes while Tj- is 
a simple convolution over the twist-2 unpolarised SPDF. T^ are more involved, and arise 
from twist-3 SDPFs, but which themselves can be expressed as convolutions of coefficient 
functions with the twist-2 SPDF. 

^We use the convention a^^^fo''^ = a^fe" + a'^fc'' and al'^fo'"! = a^'b'' ~ a'^b''. 
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To begin, we write down the tree-level OPE for the product of currents explicitly, to 
leading twist and for one flavour of quark: 

T{ny)f{z)) 



27r2 {y-zY 
i {y - z)a 
27r2 {y-zY 



^p{y)rrrHz) - Hz)rrrHy) 



X 



even j=2 w / ' 



odd 3=1 \J 



where 



(40) 



(41) 



Despite the expansion into local operators in the second step, current conservation is still 
satisfied up to equation of motion pieces. Explicitly, 

d,T{f{y)f{z)) ^ 
i (y - Z)a 



27r2 {y-zY 



odd j=l (-^ ~ ■'■)' 



- y)ni ■■■{z- y);u,._i 



even j =2 U 

- ie'^"^^ ^7^7^ i a '^^ • • • a '^^-^ ^i^ihiy + ^)) 



(42) 



= 0. 



The expression vanishes because the pairs of terms in parentheses cancel, due to the 
following identities, which are satisfied up to an equation of motion, and neglecting quark 
masses and the strong coupling: 

^ ^[M a H '^^ • • • 9 = '-e""''^ a„ i^^pYd^'---d V', (43) 



(44) 
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At first sight, it seems that it will not be possible to arrange for the quark operators 
of each parity to appear in two separate Lorentz structures, each satisfying current con- 
servation. This is not the case, however [0. The quark operators, as they appear in the 



OPE in (^) can be decomposed in terms of purely twist-2 (totally symmetric) operators, 
in a series with increasing numbers of total derivatives. To show this, we start by taking 
the quark operators and decomposing them into two pieces, of twist-2 and twist-3 (one 
pair of indices antisymmetrised) terms as follows: 



O 



o 



+ 



2(j - 1; 



S A O] 

a{ti} "/^i 



(a)Mi---Mj- 



(45) 



where the S and A operators symmetrise and antisymmetrise, respectively, the indices 
listed in subscript. Also, Oj G {O, O5} and 



Ml 



d 



s o 

a{fj.} 



(46) 
(47) 
(48) 



In order to express ( ^5] ) purely in terms of twist-2 pieces, one performs this iterative 
procedure: Substitute the twist-3 terms in ( ^5]) with eqns. (^)-(|4^); decompose the 
substituted terms again using (^5l); repeat the first two steps until there are no more 



twist-3 terms. The resulting series can be summed to give 



+ 



+ 



Qafii-fij^l 

E 

even fc=0 



1 

- s 



k+l 



i-3 

E 

odd k=l 



k+l 



{j-k- 1) d'''--- d"' e^^'^'^ dfs O 



7A'fc+2"-A'i-l 
5 



(49) 



and the same result again with O and O5 swapped. All derivatives in this expression are 
total derivatives acting on the local operators. The expression differs slightly from that 
in ref. which does not contain the term with the □, but this term can be discarded 
anyway since 

■■■ ^ = -A^d^d X-- - ^ + EOM. (50) 
Thus such terms correspond to traces and are therefore higher twist. 

By writing the matrix elements of the twist-2 operators in terms of reduced matrix 
elements, one can proceed as in the previous section and derive the form of the OPE in 
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momentum space. The matrix elements in unpolarised proton states are 



V-1 



S E ^A'^i 

■t^^ even k=0 



l/^Mfc pl^k+l ...pH (0"^'^), 



P2 



Pi 



0. 



(51) 



The matrix elements of the twist-2 operator O5 vanish because it is not possible to write 
down a completely symmetric, odd-parity Lorentz structure built only out of the vectors p 
and A. Now we can find the form of the matrix elements of the quark operators including 
twist-3 contributions. It will be convenient for our purposes to write these not in terms 
of A but of A_L, and keep only those terms up to 0(Aj_). Using eqs. (^9|) and (|51|), we 
get 



where 



P2 



o 



AMI 



M 2^ 



A 



2i 



Pi 



{m} 



(52) 
(53) 



oj(0 



of(e) 



oj(0 



even fe=0 

1 0^(0 - 



J - 1 . 



^ ^ / even k=2 
d 



even k=2 



d 



(54) 
(55) 
(56) 
(57) 



even A;=0 



All that remains for us to obtain the tree-level, large-x expansion is to take the Fourier 
transform of the off-forward matrix elements of the r.h.s. of eq. (^0]), then substitute the 
matrix elements in eqs. ( |52D and (|53|). This we write below, organised so that we can read 
off the structure functions easily: 

I J dhe''^-^ {P2\T{ra^)f{-^,z))\p,) 

(1 \ °° 

+ {p^^'q"^ + 2xp^p'' + ^p^^Afj] E OfiO 

1 ' P 1 even j=2 
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+ 



1 



+ 



2^q-p 
1 

2q ■ p 



odd j=l 

oo 

. even j=2 

(58) 

Note that the structure corresponding to Tt receives contributions only from the twist-2 
pieces of the matrix elements O^. Also, vanishes as claimed, while receive twist-3 
contributions only. These we write out explicitly so that we can check that they obey the 
factorisation condition of eq. (p8|) for DVCS: 



d 



(59) 



1 . d 



J + 



even k=2 



J + 1 



d 



even k=0 \ 

Here, the Wilson coefficients are 



1 - 5kfl, 

1 



c, 



^k,Q 



J + 1 
J + 1 



3-k + l 
J(J + 1) ' 

- (1 - 4,o) 



1 



- fl - 5, 



J(J + 1) 
j - k 



j(j + i)' 

One can check that, in each case, the coefficients obey 



3' 



(60) 

(61) 
(62) 
(63) 
(64) 

(65) 



in the large j limit. Thus there is no problem with factorisation. 



Finally, we reconstruct the factorised forms for the structure functions Tt and in 
the same way as for the simple model considered at the end of the last section. We obtain 

1 /•! du 



TT{x,^,Q\t) 



(66) 
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TLix,C,Q',t) 
T,-{x,C,Q',t) 



= 0, 



(67) 



_ 1 du 
~ ~ 2 7-1 T 



u 



2 £-2 



X 



X 



X 



1 + - In 1 - - 



u 



u 



+ 



u 



X 



ln( 1 



u 



X 



e 



2 J-i u 



1 /•! du 



d 



(68) 



2c2 



y2 _ ^2)2 



+ 



In 1 



+ 



u 

X 

u 



+ ± In 1-^) + 
uj\x 



+ {x ^ —x), 



U 



ln(l--) + 

X, 



u 



X 



U X 



V? — ^'^ u 



In 1 



u 



X 



(69) 



[X 



-x). 



In these expressions, the twist-2 quark SPDF is introduced such that 



(70) 



Notice that only its odd part w.r.t. u can contribute to any of the structure functions. 
Note also that, as claimed, the limit x ±C, is finite and that all singularities as w — > ±^ 
cancel out both for the coefficients of q{u, ^, t) and those of its derivative w.r.t. ^. 



5 Summary 

In summary, we have demonstrated factorisation for the DVCS amplitude at leading 
twist using the local OPE. We showed the three main steps in deriving this factorisation: 
(1) Obtain the large x with the OPE and identify the coefficients with the moments of 
the spectral function; (2) Invert the moments; (3) Reconstitute the factorised amplitude 
using the dispersion relation. These steps we illustrated with a simple model obeying the 
condition under which the factorisation works for the DVCS kinematics, ie. on-shell final 



17 



state photon. We also showed that all soft singularities cancel out or are integrable, given 
that the SPDFs are continuous at the end of the ERBL region. 

We then considered the tree-level OPE in QCD, including twist-3 contributions. By 
defining the structure functions for the amplitude suitably, we obtained them in terms 
of coefficient functions and the twist-2 quark SPDF. The coefficient functions do indeed 
obey the factorisation condition. 



Appendix 

This section provides a summary of the main properties of Skewed Parton Distribution 
Functions (SPDFs). We first define a twist-2 SPDF f{u,C,,t) and briefiy note some of its 
properties. Then we consider its continuity at the edges of the ERBL region — 1^| < u < 

In a scalar theory, the (twist-2) SPDF is defined as the Fourier transform of the off- 
forward matrix element of the time-ordered product of two scalar operators, 

1 

f{u, e, t) = — n / dX e*"" {p2 I T( 0(-|An) 0(iAn) ) | pi) , (71) 



27r 



oo 



where n'^ is such that = 0, n ■ p = 1 and n.A = 2^. In the limit ^ — 0, this reduces to 
the ordinary parton distribution, f{u). This limit can be taken by p2 — > Pi- 

Important properties of f{u,C,,t) follows: 

1. f{u, ^, t) = —f{—u, C,, t). This can be seen by replacing u —>■ —u in (fflD and changing 
integration variables A ^ — A to compensate. 

2. f{u,^,t) = f{u, —^,t). This follows from crossing symmetry, whereby we swap pi 
and p2. 

3. f{u,^,t) has support — 1 < n < 1. To see this, we take moments of (ITT]) as follows. 



du u^ ^ /(M,^,t) 

1 roo poo ( r/ \ 

2^ i-J"" \-J^ ' H'^X^) ) I Pi) [-'^1 e^'" 

: 0(0) id^'id''' 0(0) : 



■ ■ ■ Uf,^ [ P2 



Pi ■ (72) 
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Thus, the moments are equal to matrix elements of local, composite operators. The 
support properties of / follow from the analyticity of the these matrix elements. In 
general, if we have an analytic function g(j), and 



duu^'^ f{u) = g{j), (73) 
and provided g(j) has no poles for Rej > c, with c a positive real number, then 



fiu) = e 1- 



u 



with / a function or distribution. For a = 1, we get a function with support — 1 < m < 1. 
The matrix elements can be expected to have these properties, so that f{u,C,,t) does 
indeed have support —1 < u < 1. 

The moments of / have another important property: They are even polynomials in ^ 
of order j. To see this we decompose the matrix elements into a series of all the possible 
Lorentz structures, given that the matrix elements are those of a scalar target: 



= E C (O™-^) (75) 

even m=0 

Thus the moments are indeed polynomials of ^ of order j. 

Let us explore now the analyticity of these polynomials. If we assume that the reduced 
matrix elements {O"^'^) can be treated as analytic functions of j, then all that remains is 
to rearrange the sum over powers of ^ by expressing it as a difference of two infinite series: 

/I oo oo 

duu^-' f{u,^,t) = E r{Om,j) - e E rio""^")- (re) 
^ even m=0 even m=2 

By assumption, the first infinite sum is an analytic function of j. However, the second 
sum has a factor of . By eq. (ff^), it then follows that 



f'iT^U)- (77) 



/(n,e,t) = eii-\u\) fiu,^,t) - e\^i- 

Here, / and /' are functions or distributions of u, such that 

„1 _ oo 

/ duu^-' f{u,^,t) = E r(0'"'^), (78) 

•'• even m=0 
„! _ oo 

J duu^-' f'iu,^,t) = E r(0'"+^'^)- (79) 

^ even m=2 
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Eq. (^) is important because it explicitly separates two contributions to the SPDFs. 
The first term on the r.h.s. we shall refer to as the DGLAP contribution since its evolution 
is given by a generalised DGLAP equation. The second term, which vanishes in the 
forward limit ^ — > 0, shall be called the ERBL contribution, since it resembles, in its 
support properties, light-cone wavefuntions which have an ERBL-type evolution. The 
region |n| < |,^| region we refer to as the ERBL region. 

Eq. ( ffTl) enables us to see under what condition the SPDFs are continuous at the 
edges of the ERBL region \u\ = \^\. This is significant because integrability of the soft 
singularities which appear in the DVCS coefficient functions requires it. If f'{u/C,,^,t) 
vanishes as \u\ —>■ |^| then the SPDF is indeed continuous. A sufficient condition for this 
to happen is that the moments of /' vanish as fast as for large j. Explicitly, 

oo 1 

E r(0'"+^'^) ^ ^ f'iu/^,^,t) '"M^'o. (80) 

even m=2 3 

To our knowledge, no proof of this behaviour as yet exists, although it is believed to be 
true. 
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